Nonlinear Photoluminescence Excitation Spectroscopy of Carbon Nanotubes: 
Exploring the Upper Density Limit of One-Dimensional Excitons 
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We have studied emission properties of high-density excitons in single-walled carbon nanotubes 
through nonlinear photoluminescence excitation spectroscopy. As the excitation intensity was in- 
creased, all emission peaks arising from different chiralities showed clear saturation in intensity. Each 
peak exhibited a saturation value that was independent of the excitation wavelength, indicating that 
there is an upper limit on the exciton density for each nanotube species. We developed a theoretical 
model based on exciton diffusion and exciton-exciton annihilation that successfully reproduced the 
saturation behavior, allowing us to estimate exciton densities. These estimated densities were found 
to be still substantially smaller than the expected Mott density even in the saturation regime, in 
contrast to conventional semiconductor quantum wires. 
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Linear and nonlinear optical processes in one- 
dimensional (1-D) semiconductors arc affected by strong 
Coulomb interactions and have been the subject of many 
theoretical [J 0, 0] and experimental % 0, 0, 0, U 
Tol . 11 1 studies. Early studies of lasing in semiconductor 



quantum wires (QWRs) [B|, 0] stimulated much interest 
in the fundamental properties of high-density 1-D exci- 
tons, or correlated electron-hole (e-h) pairs, especially 
in their stability against biexciton formation, band-gap 
renormalization (BGR), and dissociation. Although it 
has been established that excitons in QWRs are stable up 
to very high densities (10 5 -10 6 cm -1 ), unanswered ques- 
tions still remain. As the density increases, an insulating 
exciton gas is expected to become unstable and eventu- 
ally transform into a metallic e-h plasma at the Mott 
density, where the inter-exciton distance approaches the 
exciton size. At what density gain should appear and 
whether a clear Mott transition exists in 1-D systems are 
open questions both theoretically 0, 0] and experi- 
mentally 0,0 0(13, EH. 

Single-walled carbon nanotubes (SWNTs) [l2| have 
recently emerged as novel 1-D solids with very strong 
quantum confinement. The exciton binding energies 
reported for semiconducting SWNTs are very large 
(~400 meV fli . fli p compared to typical GaAs QWRs 
(~20 meV [2, |6(). Despite much recent progress in un- 
derstanding their basic optical properties, so far only a 
limited number of studies have been performed under the 
condition of high carrier/exciton densities 12 1- I n 

particular, there have been no reports quantifying exci- 
ton densities in relation to the Mott density. 

Here we describe detailed characteristics of photolumi- 
nescence (PL) and photoluminescence excitation (PLE) 
spectra of SWNTs at high exciton densities. Although 
PL intensity rapidly saturated with increasing pump flu- 
ence, PL spectra were stable, indicating that there is 
no BGR, screening, or dissociation into an e-h plasma. 
However, we observed significant broadening and even- 



tual complete flattening of PLE spectra at high laser 
intensities. We show that this originates from strong 
exciton-exciton annihilation (EEA) that provides an up- 
per limit for the exciton population in the lowest-energy 
(En) state. Our model, taking into account the diffusion 
length of excitons in the EEA process, allowed us to esti- 
mate exciton densities from experimental PL saturation 
curves. The estimated densities for the highly saturated 
regime were on the order of 1 x 10 5 cm -1 , which remains 
more than one order of magnitude smaller than the ex- 
pected Mott density and explains the observed stability 
of PL spectra even in the saturation regime. 

The sample we studied was a centrifuged supernatant 
of CoMoCAT SWNTs dispersed by 1 wt% sodium cholate 
in D2O. The solution was held in a 1-mm-thick quartz cu- 
vette. The optical density in the E22 region was below 
0.2, which helped avoid inhomogencous excitation and 
re-absorption of PL. Our excitation source was a 1 kHz, 
~250 fs optical parametric amplifier (OPA), tunable in 
the visible and near-infrared, pumped by a chirped pulse 
amplifier (Clark-MXR, CPA2010). The OPA beam was 
focused onto the sample with a spot size of 300-400 /im. 
The PL from the sample was recorded with a liquid- 
nitrogen-cooled InGaAs array detector. 

Figure 1(a) compares two PL spectra. The black curve 
represents excitation by the OPA with a wavelength of 
654 nm (or 1.90 eV) and a pulse energy of 29 nJ, while 
the red curve represents excitation by a weak (100 /iW) 
CW laser with a wavelength of 658 nm (or 1.88 eV). It 
is seen that the relative intensities of different PL peaks 
are drastically different between the two curves. The in- 
set confirms that the two spectra actually coincide very 
accurately when the OPA pulse energy is kept very low 
(300 pJ). FigureQJb) shows PL spectra for pulse energies 
of 1 nJ (curves #1 and #7), 4 nJ (#2 and #6), 10 nJ (#3 
and #5), and 30 nJ (#4). The (7,5) peak is dominant at 
low fluences while the (6,5) peak becomes dominant at 
high fluences. It is important to note that the different 
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FIG. 1: (color online). Pump-intensity-dependent PL spec- 
tra, (a) Black — spectrum obtained with OPA (654 nm, 
29 nj). Red curve — spectrum obtained with a CW laser 
diode (658 nm, 100 /jW). Inset shows that the two spectra 
coincide when the OPA pulse energy is very low (300 pJ). 
(b) Change of PL spectra versus OPA pulse energy between 
1 nj and 30 nj (in the order of #1 to #7). Curve #4 corre- 
sponds to the highest fluence (~1.3 x 10 14 photons/cm 2 ). 



curves were taken in the order of #1 to #7, demonstrat- 
ing that the observed changes are reproducible and are 
not caused by any laser-induced permanent change in 
the sample. We also reproduced a similar power depen- 
dence in a dried film, showing that the power dependent 
changes are not an artifact caused by the fluidic nature 
of the sample. Additionally, note that the PL intensities 
tend to saturate at high laser fluences, while their peak 
positions do not change at all and linewidths increase 
only slightly (by ~15% from curve #1 to #4). 

Figures [2(a)-|2(c) plot integrated PL intensities for 
(6,5), (7,5), and (8,3) tubes, respectively, versus pump 
fluence up to ^1.02 x 10 14 photons/cm 2 with pump 
wavelengths of 570, 615, and 654 nm, showing clear 
saturation behaviors at high pump fluences. The E22- 
resonant wavelengths of these tube types are 570, 647, 
and 673 nm, respectively. The integrated PL intensi- 
ties were obtained through spectral decomposition anal- 
ysis by assuming 50% Lorentzian + 50% Gaussian, while 
keeping all the peak positions fixed. Figure [2] indicates 
that saturation starts at a lower (higher) fluence when the 
sample is resonantly (non-resonantly) excited. The unex- 
pectedly fast saturation of the (7,5) peak with 570 nm ex- 



FIG. 2: (color online). Integrated PL intensity versus pump 
fluence for (6,5), (7,5), and (8,3) SWNTs. Pump wavelengths 
were 570 nm (circles), 615 nm (squares), and 654 nm (trian- 
gles). The error bars account for ±5%. The solid and dashed 
curves are fitting by Eq. (1) and Eq. (2), respectively. 



citation (which is non-resonant) is likely due to its prox- 
imity to the phonon sideband at 585 nm [l8j ]. The solid 
and dashed curves are theoretical fits to the experimental 
curves based on our model to be discussed later. 

Figures E{a)-[2Id) show PLE maps taken with vari- 
ous pump fluences. The step size for the pump pho- 
ton energy was 20 meV. The data taken with the low- 
est fluence (1.2 x 10 12 photons/cm 2 ) [Fig. |3(a)] is es- 
sentially the same as that taken with low-intensity CW 
light. However, as the fluence is increased [Figs. 13(b)- 
[3(d)], the E22 excitation peaks gradually broaden and 
eventually become completely flat at the highest fluence 
(1.2 x 10 14 photons/cm 2 ) — i.e., PL intensities become 
independent of the excitation wavelength. The corre- 
sponding PLE spectra are shown in Figs. E3e)-[3(h) for 
three PL wavelengths at 983, 1034, and 1125 nm. 

To provide further insight into the nature of the ob- 
served broadening and flattening of PLE spectra, we 
also performed absorption measurements in the E22 
region using OPA pulses. Figure [J] compares two 
transmission spectra measured with the OPA at flu- 
ences of 1.0 x 10 12 photons/cm 2 (linear regime) and 
1.0 x 10 14 photons/cm 2 (saturation regime), respec- 
tively, and a transmission spectrum taken with weak CW 
whitelight. It is clear from the figure that the E22 absorp- 
tion peaks do not exhibit any broadening and shifts even 
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FIG. 3: (color online). Evolution of PLE data with increasing 
pump pulse fluence: (a) 1.2 x 10 12 , (b) 1.2 x 10 13 , (c) 4.1 x 
10 13 , and (d) 1.2 x 10 14 photons/cm 2 , (e-h): PLE spectra 
corresponding to (a)-(d) at emission wavelengths of 983 nm 
(circles), 1034 nm (squares), and 1125 nm (triangles). 



in the saturation regime, strongly suggesting that these 
intensities are not high enough to cause light-induced 
nonlinear effects such as state filling, carrier-density- 
dependent dephasing, and non-perturbative light-matter 
coupling (or "dressing") on states in the E22 range. 

We interpret these observations as results of very ef- 
ficient EE A 17, [l^|. We assume that the formation 
of En excitons occurs in a very short time scale (e.g., 
^40 fs shorter than our pulse width (^250 fs). 

Thus, excitons quickly accumulate in the En state dur- 
ing and right after photo-creation of e-h pairs. However, 
the number of excitons that can be accommodated in the 
En state is limited by EEA. As the exciton density, n x , 
approaches its maximum value, EEA begins to prevent 
a further increase by efficiently removing excitons non- 
radiatively, which explains the PL saturation behavior. 
Since EEA serves as a bottleneck for the exciton density, 
the PL intensity becomes insensitive to whether the ex- 
citons were created resonantly or non-resonantly around 
the E22 level and independent of the pump wavelength. 
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FIG. 4: (color online). Transmission spectra mea- 
sured with high-fluence OPA pulses (filled circles: 1.0 x 
10 12 photons/cm 2 , open circles: 1.0 x 10 14 photons/cm 2 ) 
compared with a transmission spectrum taken with a weak 
CW white-light beam (dashed line). 



Namely, at very high pump fluences, the PL intensity 
is determined not by how efficiently excitons are cre- 
ated but by how many En excitons can be accommo- 
dated within a particular type of SWNT as well as by the 
relative abundance of that type of SWNT in the sample. 
To determine exciton densities n x from the PL satura- 



tion curves, we have developed a theoretical model [21 1. 
Existing EEA models 171 LLSl assume the annihilaiton 
rate to be probability of finding two excitons at 

the same position in space. However, its validity becomes 
questionable when one deals with high density excitons 
in ID where the dynamics and/or size of the excitons are 
expected to influence the annihilation rate. We introduce 
a dimcnsionlcss exciton population (0 < ( < 1) defined 
as C = NI x /Lnt, where N is the number of excitons in 
a SWNT, L NT is the length of the SWNT, and l x is the 
characteristic length scale each exciton "occupies" in the 
SWNT. In a static limit, l x should simply be the exciton 
size. However, recent experiments have indicated the im- 
portance of exciton diffusion [H [H 0]. In particular, 
Cognet et al. [24| , via micro-PL studies on single micelle- 
suspended SWNTs of various chiralities, found that exci- 
tons diffusively traverse a substantial distance (~90 nm) 
within their lifetime. Therefore, we assume that l x is de- 
termined by the diffusion length and any two excitons 
formed within l x undergo EEA. We also assume that an 
exciton promoted to a higher energy level as a result of 
EEA returns to the En level with 100% probability due 
to ultrafast and efficient E22~to-En relaxation; i.e., two 
excitons become one exciton through EEA. 

Two limiting cases are considered, which we refer to as 
the "instantaneous" and "steady-state" limits. In the in- 
stantaneous limit, the initial creation of excitons is com- 
pleted before EEA follows, while the steady-state limit 
corresponds to CW excitation. The actual experimental 
situation is considered to be closer to the former. Assum- 
ing Lnt 3> ^x, the PL saturation equation is given 21 1 
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for the instantaneous limit 
and for the steady-state limit 
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These are implicit equations relating the PL intensity 
(Ipl) and pump intensity (I pump ) in terms of their re- 
spective dimensionless parameters, £ and tp. They con- 
tain no fitting parameters other than the two linear scal- 
ing factors ci and C2 and simply become a linear rela- 
tionship (ip — C) in the low density limit (£ — > 0). The 
solid and dashed curves in Fig. [5] are fits to the data by 
Eqs. (1) and (2), respectively, showing good agreement. 
The values of £ were obtained through this analysis for 
respective data points shown in Fig. [2] Finally, the exci- 
ton density can be obtained through n x = C/l x . 

Using l x = 45 nm (one half of the excursion range 
in [2J]), we estimated n x for (6,5) to be 1.7 x 10 5 , 1.3 
x 10 5 , and 1.1 x 10 5 cm" 1 for excitation with 570, 
615, and 654 nm, respectively, at a fluence of 1.02 x 
10 14 photons/cm 2 in the instantaneous limit. These 
are much smaller than the expected Mott density n* 
(~7 x 10 6 cm -1 , assuming exciton size ~ 1.5 nm [25|), 
which explains the stability of PL (Fig. [Tj in the satu- 
rated regime. However, the densities of as-created e-h 
pairs is estimated to be 1-2 x 10 6 cm -1 in the cases of 
resonant excitation, which is similar to n*. These obser- 
vations appear qualitatively different from GaAs QWRs 
where PL saturation is not obvious until the formation 
of biexcitons and an e-h plasma (ll| . Highly efficient and 
rapid EE A in SWNTs, which is consistent with the re- 
ported absence of biexciton signatures [26| , is probably 
the direct reason for such characteristic differences. 

In summary, we have studied PL and PLE spectra of 
SWNTs at high exciton densities. Complete flattening 
of PLE spectra and clear PL saturation were observed, 
indicating the existence of an upper density limit. We 
developed a model that reproduced the observed satura- 
tion behavior and allowed us to estimate exciton densi- 
ties, which remained more than one-order-of-magnitude 
lower than the expected Mott density and explained the 
stability of PL spectra in the saturated regime. 
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